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Electron  radiation  induced  hydrocarbon  contamination  can  be  either  a  problem  or  a  useful  tool  in 
electron  beam  analyses  and  lithographies.  We  have  used  atomic  force  microscopy  to  study  electron 
beam  written  contamination  structures.  Contamination  is  shown  not  only  to  arise  from  the  primary 
electron  beam  but  also  from  the  energy  scattered  outside  of  the  direct  impingement  area.  The  size 
of  the  contamination  structures  correlates  well  with  that  expected  from  electron  scattering  theory.  By 
varying  the  geometry  of  the  written  structures,  the  rate  at  which  the  electron  dose  is  deposited,  and 
the  nearby  hydrocarbon  surface  density,  we  show  that  surface  diffusion  of  hydrocarbon  molecules 
plays  a  primary  role  in  the  formation  of  the  contamination  structures.  ©  1996  American  Vacuum 
Society. 


I.  INTRODUCTION 

It  is  well  known  that  exposure  of  a  surface  in  vacuum  to 
an  electron  beam  can  lead  to  the  buildup  of  a  carbon-rich 
film  in  the  exposed  area.1  This  contamination  is  produced 
through  electron  beam  induced  polymerization  of  hydrocar¬ 
bon  molecules  that  are  adsorbed  onto  the  surface.  The 
sources  of  these  hydrocarbon  molecules  include  residual  hy¬ 
drocarbons  in  the  vacuum  chamber  and  layers  adsorbed  onto 
the  surface  during  sample  preparation.2  The  contamination 
can  degrade  specimen  images  in  transmission  and  scanning 
electron  microscopies.3  Care  must  also  be  taken  to  avoid  this 
contamination  in  certain  electron  beam  lithographies.  For  ex¬ 
ample,  patterning  of  oxidized  silicon  has  been  achieved  by 
directly  writing  on  the  substrate  with  an  electron  beam  and 
then  partially  etching  the  oxide.4  The  presence  of  a  contami¬ 
nation  layer  in  the  exposed  region  can  interfere  with  this  etch 
process.  Sacrificial  layers,5  and  pre-  and  postexposure  sub¬ 
strate  cleaning6-8  have  been  used  to  address  this  problem. 
Other  electron  beam  lithographies  where  contamination  is  an 
issue  include  selective  etching  of  GaAs9  and  selective  area 
growth.10  However,  contamination  is  not  always  undesirable. 
Electron  beam  written  contamination  structures  have  been 
used  as  masks  for  pattern  transfer. 1 1  Metal  lines  with  widths 
down  to  80  A  have  been  fabricated  using  contamination 
masks.12  Contamination  is  also  utilized  in  film  thickness 
measurements  made  in  transmission  electron  microscopy.13 

Atomic  force  microscopy  (AFM)  is  a  valuable  technique 
for  analyzing  the  topography  of  small  structures.  In  this  ar¬ 
ticle  we  present  the  results  of  an  AFM  study  of  electron 
beam  written  hydrocarbon  contamination  structures.  We  dis¬ 
cuss  the  size  and  shape  of  a  typical  structure  and  relate  the 
lateral  extent  of  the  contamination  to  electron  scattering 
theory.  The  dependence  of  the  developed  contamination  on 
substrate  material,  electron  dose,  delay  time  between  mul¬ 
tiple  exposures  (average  electron  deposition  rate),  and  writ- 
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ten  geometry  is  presented.  Finally,  a  surface  diffusion  picture 
is  developed  that  models  the  observed  dependencies  and  pro¬ 
vides  an  estimate  for  the  characteristic  time  of  the  hydrocar¬ 
bon  diffusion. 

II.  EXPERIMENTAL  PROCEDURE 

The  hydrocarbon  contamination  structures  were  written  in 
a  JEOL  6400  scanning  electron  microscope  modified  for  li¬ 
thography.  The  exposures  were  made  at  an  accelerating  volt¬ 
age  of  40  kV  and  a  beam  current  of  100  nA.  Electron  doses 
were  controlled  by  the  duration  of  an  individual  exposure 
and  by  the  number  of  times  the  region  was  exposed.  The 
specimen  chamber  was  maintained  at  a  pressure  of  2X10  7 
Torr  dining  the  exposures. 

Various  substrates  were  used  in  this  study.  These  include 
semi-insulating  and  n+  (100)  GaAs,  (100)  GaAs  precoated 
with  100  nm  of  Al06Ga04As  through  metalorganic  chemical 
vapor  deposition,  and  p -doped  (100)  Si.  The 
Al0  6Ga0  4As/GaAs  samples  were  grown  for  a  separate  study 
of  InAs  selective  area  growth  on  GaAs  in  which  the 
Al06Ga04As  layer  would  be  later  etched  away  prior  to  the 
InAs  growth  in  order  to  remove  the  electron  beam  induced 
contamination.10 

After  the  electron  beam  exposure,  the  contamination 
structures  were  analyzed  under  ambient  conditions  using  a 
Park  Scientific  Instruments  atomic  force  microscope  (AFM). 
Geometrical  information  was  obtained  from  contact  mode 
AFM  images.  The  lateral  sizes  of  the  structures  were  typi¬ 
cally  of  the  order  of  10  pm,  while  the  heights  were  at  most 
200  A.  Measuring  such  a  small  height  variation  over  a  large 
lateral  distance  required  a  nonlinear  correction  during  image 
acquisition  to  compensate  for  the  AFM  scanner  nonlineari¬ 
ties.  Additionally,  for  the  images  in  which  a  significant  un¬ 
modified  region  of  the  substrate  was  present,  a  second  cor¬ 
rection  was  performed.  The  correction  consisted  of  the 
subtraction  from  the  entire  image  of  a  third-order  polyno¬ 
mial,  fit  to  the  unmodified  and  presumably  flat  regions  of  the 
substrate.  No  other  image  modifications  were  made. 
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Fig.  1.  (a)  Atomic  force  microscope  image  of  an  electron  beam  written 
contamination  line.  The  written  line  is  1  mm  in  length  of  which  a  section 
near  the  middle  of  the  line  is  shown  in  this  image.  The  contamination  struc¬ 
ture  was  created  by  passing  the  electron  beam  over  the  same  vertical  line 
(labeled  1)  40  times  for  a  total  line  dose  of  367  nC//xm.  The  regions  modi¬ 
fied  indirectly  by  the  electron  beam  and  unchanged  by  the  electron  beam  are 
labeled  2  and  3,  respectively.  The  substrate  is  n+  (100)  GaAs.  The  scan  size 
is  30  yumX30  //m  with  a  total  height  variation  of  151  A.  The  dashed  line 
indicates  where  the  cross-sectional  profile  shown  in  (b)  was  taken,  (b) 
Cross-sectional  profile  of  the  electron  beam  written  line  shown  in  (a).  Be¬ 
tween  the  dashed  lines  is  the  region  in  which  the  scattered  electrons  interact 
with  the  substrate  surface  as  determined  from  a  modified  electron  diffusion 
model. 


III.  RESULTS  AND  DISCUSSION 

An  AFM  image  of  a  typical  contamination  structure  writ¬ 
ten  with  the  electron  beam  is  shown  in  Fig.  1(a).  The  image 
is  of  a  30  fjm  section  of  a  1  mm  long  vertical  line  on  a  GaAs 
substrate.  Three  distinct  regions  are  evident  from  Fig.  1:  (1) 
a  relatively  sharp  center  peak,  (2)  a  broader  surrounding 
raised  area,  and  (3)  the  unmodified  substrate.  The  first  region 
is  in  the  area  of  direct  electron  beam  impingement.  This 
primary  beam  of  electrons  provides  the  energy  necessary  to 
polymerize  the  normally  mobile  surface  hydrocarbons,  cre¬ 
ating  an  immobile  hydrocarbon  contamination  structure  the 
width  of  the  electron  beam.  However,  the  interaction  of  the 
primary  beam  alone  does  not  fully  explain  the  observed 
structure.  As  the  primary  electrons  enter  the  substrate,  a 
number  of  interactions  take  place.14  One  such  interaction  is 
the  scattering  of  the  primary  electrons  by  the  atomic  elec¬ 
trons  and  nuclei  of  the  substrate,  which  cause  the  electron 
beam-sample  interaction  region  width  to  be  larger  than  the 
beam  diameter.  In  addition,  the  primary  electrons  may  col¬ 
lide  with  the  atomic  electrons  and  ionize  them,  creating  sec¬ 
ondary  electrons  that  will  subsequently  scatter  about  in  the 
substrate.  Some  of  the  scattered  primary  and  secondary  elec¬ 
trons  pass  through  the  substrate  surface  causing  further  hy¬ 
drocarbon  polymerization.  These  interactions  give  rise  to  the 
broad  contamination  feature  of  region  2.  As  might  be  ex¬ 
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Fig.  2.  Schematic  cross  section  of  an  electron  beam  exposed  solid.  The 
primary  electrons  that  have  entered  the  solid  are  scattered  by  the  atomic 
electrons  and  nuclei  of  the  solid.  The  primary  electrons  can  also  ionize  the 
atomic  electrons  giving  rise  to  secondary  electrons  that  scatter  about  in  the 
solid.  The  interaction  sphere  is  the  approximate  volume  in  which  the  pri¬ 
mary  and  scattered  electrons  have  energies  greater  than  their  thermal  ener¬ 
gies.  The  depth  of  this  sphere  into  the  solid  is  the  maximum  electron  range 
R,  and  the  radius  of  the  intersection  circle  between  the  interaction  sphere 
and  the  surface  of  the  solid  is  the  backscattering  radius  rB  .  Both  the  incident 
primary  electrons  and  the  electrons  scattered  back  through  the  surface  can 
lead  to  hydrocarbon  polymerization  on  the  surface  of  the  solid. 


pected,  the  contamination  is  greater  in  the  primary  beam  area 
and  decreases  with  distance  from  this  area  due  to  the  smaller 
amount  of  energy  being  deposited  away  from  the  direct  im¬ 
pingement  region.  Region  3  is  beyond  the  area  of  significant 
electron  beam-sample  interaction  and  therefore  exhibits  no 
electron  beam  induced  contamination. 

The  spatial  extent  over  which  the  scattered  electrons  in¬ 
teract  with  the  substrate  surface  is  indicated  by  the  size  of 
region  2.  This  experimentally  determined  interaction  length 
can  be  compared  to  values  obtained  from  electron  scattering 
theory.  The  electron  interaction  volume  in  the  substrate  is 
approximated  by  a  sphere  that  is  partially  contained  within 
the  substrate  as  shown  in  Fig.  2.  Within  this  volume,  the 
primary  and  scattered  electrons  have  energies  greater  than 
their  thermal  energies.  The  quantity  to  be  compared  to  the 
size  of  region  2  is  the  radius  of  the  circle  of  intersection 
between  the  interaction  sphere  and  the  substrate  surface.  This 
is  the  backscattering  range  rB .  In  the  modified  electron  dif¬ 
fusion  model  proposed  by  Kanaya  and  Okayama,15  the  back- 
scattering  range  is  given  by 


rB  = 


CRy 
1  +  y ' 


(1) 


where  C  is  a  parameter  whose  best  fit  value  was  determined 
to  be  1.1  in  Ref.  15.  The  maximum  electron  range  R  (/tm) 
and  the  parameter  y  of  Eq.  (1)  are 

2.76X  1(T7A£5/3  (1  +  0.978X  1(T6£)5/3 
R~  pZm  ( 1  +  1.957X  10~6£)4/3  ’  ® 

y=0.187Z2/3,  (3) 

where  Z,  A  (g),  and  p  (g/cm3)  are  the  atomic  number,  atomic 
mass,  and  density  of  the  substrate,  and  E  (eV)  is  the  energy 
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Fig.  3.  (a)  Atomic  force  microscope  image  of  an  electron  beam  written 
contamination  line.  The  written  line  is  1  mm  in  length  of  which  a  section 
near  the  middle  of  the  line  is  shown  in  this  image.  The  contamination  struc¬ 
ture  was  created  by  passing  the  electron  beam  over  the  same  vertical  line  40 
times  for  a  total  line  dose  of  367  nC/pm.  The  substrate  is  p -doped  (100)  Si. 
The  scan  size  is  40  //mX40  pm  with  a  total  height  variation  of  143  A.  The 
dashed  line  indicates  where  the  cross-sectional  profile  shown  in  (b)  was 
taken,  (b)  Cross-sectional  profile  of  the  electron  beam  written  line  shown  in 
(a).  Between  the  dashed  lines  is  the  region  in  which  the  scattered  electrons 
interact  with  the  substrate  surface  as  determined  from  a  modified  electron 
diffusion  model. 


of  the  incident  electrons.  To  apply  this  model  to  a  GaAs 
substrate  we  have  used  atomic  number  and  mass  values  that 
are  averages  of  the  individual  elemental  values.  Using  Eqs. 
(1)— (3)  with  Z=  32,  3  =  72.32  g,  p= 5.32  g/cm3,  and  £  =  40 
keV,  we  obtain  a  value  for  the  backscattering  range  of  5.6 
pm.  This  range  is  compared  to  the  cross-sectional  profile  of 
the  contamination  line  in  Fig.  1(b).  A  dashed  line  is  drawn  on 
each  side  of  the  central  contamination  region,  a  distance  rB 
from  the  center  peak.  These  lines  approximately  coincide 
with  the  outer  edges  of  region  2. 

The  backscattering  range  depends  in  part  on  the  physical 
properties  of  the  substrate,  as  seen  in  Eq.  (2).  To  further 
confirm  the  dependence  of  the  contamination  structure  on  the 
backscattering  range,  an  electron  beam  exposure  identical  to 
that  done  on  the  GaAs  substrate  was  made  on  a  Si  substrate. 
An  AFM  image  of  the  resulting  contamination  line  is  shown 
in  Fig.  3(a).  From  Eqs.  (l)-(3)  the  range  for  Si  is  calculated 
to  be  8.2  pm,  which  is  larger  than  that  of  GaAs  as  a  result  of 
the  smaller  density  of  Si.  Hence  a  wider  contamination  line 
is  expected  and  observed.  The  calculated  value  compares 
fairly  well  with  the  extent  of  the  contamination  as  shown  in 
Fig.  3(b).  Additionally,  the  height  of  the  contamination  out¬ 
side  the  direct  electron  beam  impingement  region  is  smaller 
than  that  observed  on  the  GaAs  sample.  This  is  consistent 
with  the  energy  of  the  electron  beam  being  scattered  over  a 
larger  volume  in  the  Si  sample.  The  energy  density  at  the 


Fig.  4.  Atomic  force  microscope  images  and  cross-sectional  profiles  of  elec¬ 
tron  beam  written  contamination  lines  written  with  various  electron  doses. 
Each  line  is  1  mm  long  of  which  a  10  /mnX30  gm  section  near  the  middle 
of  the  line  is  shown.  The  electron  dose  was  varied  by  changing  the  number 
of  passes  the  electron  beam  made  over  the  line.  The  substrate  is  n+  (100) 
GaAs.  The  number  of  electron  beam  passes,  total  line  dose,  and  total  height 
variation  for  each  image  are:  (a)  40  passes,  367  nC/ fjm,  149  A,  (b)  20 
passes,  183  nC//xm.  82  A,  (c)  10  passes,  92  nC//um,  46  A,  (d)  5  passes,  46 
nC//rm,  32  A.  (e)  Cross-sectional  profiles  of  the  electron  beam  written  lines 
with  the  top  through  bottom  profile  corresponding  to  images  (a)-(d),  respec¬ 
tively.  Each  profile  is  an  average  of  five  individual  line  sweeps. 


surface  is  less,  and  therefore  a  smaller  accumulation  of  con¬ 
tamination  results. 

To  achieve  high  electron  doses,  the  electron  beam  was 
passed  repeatedly  over  the  exposed  region.  The  dose  depen¬ 
dence  of  the  contamination  structures  may  be  conveniently 
studied  by  varying  the  number  of  electron  beam  passes  over 
a  given  line.  The  results  of  such  an  experiment  on  a  GaAs 
substrate  are  shown  in  Fig.  4.  The  images  in  Figs.  4(a)-4(d) 
are  of  1  mm  long  contamination  lines  with  the  number  of 
electron  beam  passes  successively  halved  in  going  from  40 
passes  in  Fig.  4(a)  to  5  passes  in  Fig.  4(d).  The  three  distinct 
regions  identified  previously  in  Fig.  1  are  visible  in  each  of 
the  lines  of  Fig.  4.  Moreover,  from  the  cross-sectional  pro¬ 
files  shown  in  Fig.  4(e),  it  is  evident  that  the  lateral  size  and 
overall  shape  of  the  regions  remain  relatively  unchanged  as 
the  number  of  passes  is  varied.  Each  pass  of  the  electron 
beam  simply  adds  approximately  the  same  amount  of  con¬ 
tamination  to  the  structure  as  the  previous  individual  passes, 
leading  to  the  height  of  the  contamination  being  proportional 
to  the  number  of  passes.  This  dependence  has  also  been  ob¬ 
served  for  similar  lines  written  on  Si  substrates. 

The  contamination  structures  discussed  to  this  point  have 
all  been  relatively  long,  1  mm  lines.  A  change  in  the  cross- 
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Fig.  5.  Atomic  force  microscope  images  and  cross-sectional  profiles  of  elec¬ 
tron  beam  contamination  lines  of  various  lengths.  Each  contamination  struc¬ 
ture  was  created  by  passing  the  electron  beam  over  the  same  vertical  line  40 
times  for  a  total  line  dose  of  367  n(7/ani.  A  10  /xmX25  gm  section  near  the 
middle  of  each  line  is  shown  in  the  images.  The  substrate  is  100  nm  of 
Al0  6Ga0  4As  on  semi-insulating  (100)  GaAs.  The  line  length,  time  delay  for 
the  beam  to  return  to  the  same  point  on  the  line  Td ,  and  the  total  height 
variation  for  each  image  are:  (a)  1  mm,  100.2  s,  137  A,  (b)  300  /a in.  30.2  s, 
100  A,  (c)  200  jirm,  20.2  s,  69  A,  and  (d)  50  jum,  5.2  s,  37  A.  (e)  Cross- 
sectional  profiles  of  the  electron  beam  written  lines  with  the  top  through 
bottom  profile  corresponding  to  images  (a)-(d),  respectively.  Each  profile  is 
an  average  of  five  individual  line  sweeps. 


Fig.  6.  Atomic  force  microscope  images  and  cross-sectional  profiles  of  elec¬ 
tron  beam  written  dots.  Each  dot  structure  was  created  by  exposing  the  same 
spot  40  times  for  a  total  dose  of  33.6  nC  with  a  delay  time  Td  between  each 
of  the  exposures.  An  8.5  /imX8.5  /arn  region  around  each  dot  is  shown  in 
the  images.  The  substrate  is  100  nm  of  Al06Ga,)4As  on  semi-insulating  (100) 
GaAs.  The  time  delay  between  the  exposures  and  the  total  height  variation 
for  each  image  are:  (a)  2.3  s,  89  A,  (b)  1.2  s,  46  A,  (c)  0.7  s.  22  A,  and  (d) 
22  ms,  21  A.  (e)  Cross-sectional  profiles  of  the  electron  beam  written  dots 
with  the  top  through  bottom  profile  corresponding  to  images  (a)-(d),  respec¬ 
tively. 


sectional  profile  of  the  line  occurs  when  the  length  of  the 
written  line  is  varied.  In  Fig.  5  a  set  of  contamination  lines  of 
various  lengths  written  on  a  Al0  6Ga0  4As/GaAs  substrate  is 
shown.  The  images.  Figs.  5(a)-5(d),  are  of  identically  writ¬ 
ten  lines,  except  that  the  line  lengths  are  varied  from  the 
longest  of  1  mm  in  Fig.  5(a)  to  the  shortest  of  50  jum  in  Fig. 
5(d).  In  the  images  of  the  longer  line  lengths  of  Figs.  5(a)- 
5(c),  the  direct  impingement  area,  region  1,  is  clearly  distin¬ 
guishable  from  the  surrounding  broader  contamination  of  re¬ 
gion  2.  However,  as  can  be  seen  in  the  cross-sectional 
profiles  of  Fig.  5(e),  the  height  of  the  structure  decreases  as 
the  line  is  shortened.  This  decrease  is  not  uniform,  but  is 
greatest  in  region  1  and  lessens  in  moving  outward  into  re¬ 
gion  2.  The  height  decrease  near  the  center  of  the  line  cross 
section  continues  as  the  line  is  shortened  until  there  appears 
to  be  little  or  no  hydrocarbon  contamination  in  this  region,  as 
is  the  case  for  the  50  /rm  line  of  Fig.  5(d).  Significant  con¬ 
tamination  is  still  present  in  region  2  away  from  the  center  of 
the  line  cross  section.  Similar  results  have  been  obtained  on 
GaAs  without  the  Al()  6Ga0  4As  layer. 

The  change  in  profile  with  line  length  is  a  direct  result  of 
the  method  in  which  the  lines  are  written  and,  additionally, 
provides  information  on  the  processes  involved  in  the  forma¬ 


tion  of  the  hydrocarbon  contamination.  As  stated  previously, 
the  structures  are  written  by  passing  the  electron  beam  over 
the  same  line  a  number  of  times  in  order  to  achieve  high 
electron  doses  and,  as  is  evident  from  Fig.  4,  to  create  sig¬ 
nificant  contamination.  The  electron  beam  is  swept  at  the 
same  constant  velocity  along  the  line  for  all  line  lengths; 
therefore  the  length  of  the  line  dictates  the  amount  of  time 
that  will  elapse  before  the  electron  beam  returns  to  the  same 
spot  on  the  line.  This  delay  time  rd  is  large  for  a  long  line 
and  small  for  a  short  one.  Figure  5  clearly  indicates  that  the 
structure  does  not  simply  depend  on  the  overall  dose,  but 
also  on  Td ,  in  other  words,  on  the  rate  at  which  the  electron 
beam  dose  is  deposited. 

Another,  simpler  structure  used  to  study  the  effect  of  the 
delay  time  is  an  electron  beam  written  dot.  A  dot  in  this  case 
is  written  by  electron  beam  exposing  a  single  spot  on  the 
substrate  for  a  short  time  (approximately  8  ms),  waiting  a 
variable  length  of  time  Td ,  and  then  reexposing  the  same 
spot.  This  process  is  repeated  until  the  desired  electron  dose 
is  achieved.  Shown  in  Fig.  6  is  a  set  of  dots,  each  written  by 
exposing  a  spot  on  an  Al0  6Ga0  4As/GaAs  substrate  40  times 
with  rd  being  successively  decreased  from  2.3  s  for  the  dot 
shown  in  Fig.  6(a)  to  22  ms  for  the  dot  shown  in  Fig.  6(d). 
The  long  time  delay  contamination  dot  of  Fig.  6(a)  consists 
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line  length  (mm) 


delay  time  (s) 

Fig.  7.  Contamination  height  of  the  direct  impingement  region  plotted  as  a 
function  of  delay  time  for  both  the  dot  and  the  line  structures.  Each  structure 
was  created  by  exposing  the  same  pattern  forty  times  for  a  total  dose  of  33.6 
nC  for  the  dots  and  367  nC/fim  for  the  lines.  The  substrate  is  100  nm  of 
Al06Ga04As  on  semi-insulating  (100)  GaAs.  The  solid  lines  are  fits  to  the 
data  of  a  model  based  on  hydrocarbon  surface  diffusion. 


of  a  relatively  sharp  center  peak  in  the  area  of  direct  electron 
beam  impingement  and  a  broader  surrounding  contamination 
region  resulting  from  the  scattered  energy  of  the  primary 
beam.  As  the  time  delay  is  made  smaller,  the  height  of  the 
structure  decreases,  as  seen  in  Figs.  6(b)  and  6(c).  From  the 
cross-sectional  profiles  in  Fig.  6(e),  it  is  evident  that  this 
reduction  of  height  is  greatest  at  the  dot  center,  which  is 
similar  to  the  effect  observed  with  the  line  structures. 

The  height  reduction  of  the  contamination  structures  with 
decreasing  delay  time  is  a  direct  result  of  the  finite  time 
required  for  the  exposed  region,  reduced  in  mobile  hydrocar¬ 
bon  molecules,  to  be  replenished  with  molecules.  In  other 
words,  the  electron  beam  exposure  temporarily  decreases  the 
number  of  mobile  molecules  in  the  exposed  region  by  poly¬ 
merizing  the  hydrocarbons  and  creating  the  observed  con¬ 
tamination.  Thus,  after  the  exposure,  the  concentration  of 
mobile  hydrocarbons  in  the  exposed  region  is  below  the 
equilibrium  value.  Hydrocarbons  must  then  be  transported  to 
this  region  in  order  to  reestablish  the  equilibrium  concentra¬ 
tion.  The  transport  process  requires  a  finite  amount  of  time 
dictated  by  the  mechanism  of  the  transport.  If  the  time  delay 
is  much  shorter  than  the  transport  time,  the  exposed  region 
becomes  depleted  of  hydrocarbons  and  the  structure  will  not 
grow  with  subsequent  exposures.  As  the  delay  time  is  in¬ 
creased,  however,  more  hydrocarbons  will  be  present  in  the 
region  of  exposure  prior  to  the  event,  leading  to  a  contami¬ 
nation  structure  of  greater  height.  If  the  time  between  expo¬ 
sures  is  much  larger  than  the  characteristic  transport  time,  the 
hydrocarbon  supply  will  replenish  between  exposures  and 
the  structure  will  grow  substantially  with  each  additional  ex¬ 
posure.  Increasing  the  delay  time  further  cannot  increase  the 
quantity  of  hydrocarbons  beyond  the  completely  replenished 
equilibrium  value.  The  contamination  height,  being  propor¬ 
tional  to  the  quantity  of  hydrocarbons  present  during  the  ex¬ 
posure,  will  therefore  also  not  increase  with  delay  time  in 
this  limit.  Qualitatively  this  is  what  is  observed  in  the  line 


and  dot  structures  of  Figs.  5  and  6.  Plots  of  the  center  con¬ 
tamination  height  as  a  function  of  the  delay  time  for  dot  and 
line  structures  are  shown  in  Fig.  7.  The  structures  were  all 
written  on  the  same  Al0  6Ga0  4As/GaAs  substrate.  As  ex¬ 
pected,  the  height  of  the  structures  increases  with  the  delay 
time.  The  height  also  changes  rapidly  for  small  delay  times 
and  then  appears  to  saturate  at  large  times,  consistent  with 
the  model. 

The  dependence  of  the  contamination  height  on  the  delay 
time  can  be  used  to  extract  information  on  the  hydrocarbon 
transport  mechanism.  Two  mechanisms  by  which  the  hydro¬ 
carbons  can  be  supplied  to  the  exposed  region  are  (1)  diffu¬ 
sion  along  the  substrate  from  an  unexposed  region,  and  (2) 
adsorption  from  the  vacuum.2  If  surface  diffusion  is  the  pri¬ 
mary  mechanism,  then  the  concentration  of  mobile  hydrocar¬ 
bon  molecules  can  be  described  by  the  following  diffusion 
equation: 


dn{x,t ) 

Jt 


=  DV2n(x,t), 


(4) 


where  n(x,t )  is  the  surface  density  of  the  mobile  hydrocar¬ 
bon  molecules  at  time  t  and  location  x,  and  D  is  the  diffusion 
constant.  The  hydrocarbon  density  dictates  in  part  the 
amount  of  contamination  that  develops  during  the  exposure. 
Given  the  density  distribution  after  a  particular  exposure 
n(x,  t  =  0),  the  distribution  that  will  be  present  at  the  start  of 
the  subsequent  exposure  n(x,  t=rd )  can  be  obtained  from 
Eq.  (4).  For  simplicity  consider  the  situation  in  which  a  re¬ 
gion  of  the  substrate  is  reduced  in  hydrocarbons  from  the 
equilibrium  density  nQ  to  the  value  nr  as  a  result  of  the 
exposure.  Applying  this  to  both  the  dot  and  the  line  geom- 
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where  2d  is  the  effective  size  of  the  exposed  region,16 

and 

polar  coordinates  x  =  (r,9)  and  rectangular  coordinates 
x  =  (x,y)  have  been  used  for  the  dot  geometry  and  the  line 
geometry,  respectively.  The  solutions  to  the  diffusion  equa¬ 
tion  (4)  assuming  these  initial  conditions  are  given  in  the 
Appendix.  In  order  to  compare  the  two  geometries,  consider 
the  time  dependence  of  the  density  at  the  center  of  the  ex¬ 
posed  region,  where  the  minimum  density  occurs.  Letting 
t=  rd ,  the  following  are  obtained  from  Eqs.  (A4)  and  (A10): 


0,  t=rd)=nr-\-(riQ— nr)exp^ 

4  DTdl 

J  dot, 

(7) 

0,  t=r(i)=nrJr{nQ— /2r)erfcj 

d  ) 

line. 

(8) 

The  two  solutions  are  plotted  as  a  function  of  the  dimension¬ 
less  delay  time  in  Fig.  8  for  a  region  that  is  initially  com¬ 
pletely  depleted,  nr=  0.  These  results,  based  solely  on  a  sur¬ 
face  diffusion  source,  predict  the  delay  time  dependency 
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Fig.  8.  Normalized  hydrocarbon  surface  density  plotted  as  a  function  of  the 
dimensionless  delay  time  for  the  dot  and  the  line  geometries.  The  normal¬ 
ized  density  at  the  center  of  an  initially  depleted  region  of  size  2d  is  plotted. 
The  change  in  density  has  been  determined  assuming  only  diffusion  of  hy¬ 
drocarbons  along  the  surface.  Both  normalized  densities  approach  unity  as 
the  delay  time  gets  large,  however  this  occurs  more  rapidly  in  the  dot  ge¬ 
ometry  due  to  the  larger  hydrocarbon  source  available  to  the  dot. 


observed  in  the  data  of  Fig.  7.  At  small  delay  times  for  both 
structures,  Td<d2/4D ,  the  density  is  small  and  therefore  the 
accumulated  contamination  is  small.  For  large  delay  times, 
rd>d2/4D ,  the  density  is  near  the  equilibrium  value,  allow¬ 
ing  significant  contamination  to  develop  during  the  exposure. 
For  intermediate  delay  times,  the  measured  height  changes 
more  rapidly  for  the  dot  geometry.  This  effect  is  also  pre¬ 
dicted  by  the  simple  surface  diffusion  model  as  the  plots  of 
Fig.  8  illustrate.  The  difference  between  the  two  geometries 
results  from  the  dot  effectively  having  a  larger  source  of 
hydrocarbons  to  draw  from  and  therefore  recovering  more 
quickly  from  the  depleted  state.  This  geometrical  effect  can 
be  further  illustrated  by  looking  at  the  end  of  the  65  /jl m  long 
line,  shown  in  Fig.  9(a).  The  end  of  the  line,  much  like  the 
dot,  effectively  has  a  greater  source  of  hydrocarbons  than  a 
segment  far  from  the  ends.  Consequently,  the  line  end  has  a 
greater  contamination  height  as  can  be  seen  from  the  cross- 
sectional  profiles  of  Fig.  9(b).  The  delay  time  dependence  of 
the  line  height  away  from  the  line  ends  is  compared  to  that  at 
the  line  ends  in  Fig.  10.  The  greater  height  rise  with  delay 
time  observed  for  the  line  end  is  consistent  with  the  geo¬ 
metrical  difference  between  the  two  regions.  These  geometri¬ 
cal  differences  would  not  be  present  if  the  hydrocarbons 
were  supplied  primarily  by  direct  adsorption  from  the 
vacuum.  In  this  simple  picture,  such  a  process  would  evenly 
coat  the  exposed  regions. 

From  the  simple  diffusion  model,  we  can  get  an  estimate 
for  the  diffusion  time  parameter  d2/4D  by  fitting  Eqs.  (7)  and 
(8)  to  the  data  of  Fig.  7.  The  best-fit  curves,  assuming  nr=  0, 
are  shown  in  Fig.  7.  The  value  of  d2/4D  is  3  s  for  the  dots 
and  10  s  for  the  lines.  The  difference  between  these  two 
values  can  be  partially  accounted  for  by  a  difference  in  the 
size  of  the  depleted  hydrocarbon  region.19  The  line  exposures 
have  been  effectively  written  with  a  greater  dose,  due  to  the 
method  of  exposure,  and  hence  the  scattered  electron  beam 
region  has  more  contamination  than  that  seen  in  the  dot 


Fig.  9.  (a)  Atomic  force  microscope  image  of  an  electron  beam  written 
contamination  line.  The  written  line  is  65  /zm  in  length  of  which  a  section  at 
the  end  of  the  line  is  shown  in  this  image.  The  contamination  structure  was 
created  by  passing  the  electron  beam  over  the  same  nearly  vertical  line  40 
times  for  a  total  line  dose  of  472  nC//xm.  The  substrate  is  n+  (100)  GaAs. 
The  scan  size  is  25  /xmX25  /zm  with  a  total  height  variation  of  74  A.  The 
dashed  lines  indicate  where  the  cross-sectional  profiles  shown  in  (b)  were 
taken,  (b)  Cross-sectional  profiles  of  the  electron  beam  written  line  shown  in 
(a).  The  cross  sections  show  that  there  is  a  greater  quantity  of  contamination 
at  the  line  ends  as  compared  to  regions  away  from  the  ends. 


structures.  The  value  of  2 d  is  then  greater  for  the  line  struc¬ 
tures. 

If  the  contamination  growth  is  governed  by  the  surface 
diffusion  process  as  suggested,  it  should  be  possible  to 
modify  the  growth  by  changing  the  hydrocarbon  density  in 


line  length  (pm) 


Fig.  10.  Contamination  height  of  the  direct  impingement  region  of  a  line 
plotted  as  a  function  of  the  delay  time  and  line  length.  The  height  at  a  line 
end  is  compared  to  that  toward  the  middle.  Each  of  the  line  structures  was 
written  by  passing  the  electron  beam  over  the  same  line  40  times  for  a  total 
line  dose  of  472  nC//zm.  The  substrate  is  n+  (100)  GaAs.  The  line  end 
height  rises  more  rapidly  with  delay  time  due  to  the  larger  source  of  hydro¬ 
carbons  available  to  the  ends. 
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Fig.  11.  (a)  Cross-sectional  profiles  of  electron  beam  written  contamination 
dots  obtained  with  an  atomic  force  microscope.  Each  dot  structure  was 
written  by  exposing  the  same  spot  40  times  for  a  total  dose  of  33.6  nC.  The 
delay  time  between  each  exposure  was  100.2  s.  From  top  to  bottom,  the 
profiles  correspond  to  an  isolated  dot,  and  dots  written  inside  concentric 
squares  with  the  smallest  square  having  a  side  of  length  80.  60,  and  40  /rm, 
respectively.  The  specific  patterns  that  were  written  are  shown  schematically 
in  (c).  The  substrate  is  n+  (100)  GaAs.  (b)  Cross-sectional  profiles  of  dots 
written  in  the  same  manner  as  those  of  (a).  The  substrate  is  /) -doped  (100) 
Si. 


the  vicinity  of  the  written  structure  prior  to  exposure.  The 
following  experiment  was  performed  to  test  this  idea.  The 
perimeters  of  a  series  of  concentric  squares  were  first  written 
in  order  to  reduce  the  hydrocarbon  density  surrounding  a 
region  on  the  substrate.  Then  a  single  dot  was  written  in  the 
center  of  the  squares.  The  entire  process  was  repeated  until 
the  desired  electron  dose  at  the  dot  was  achieved.  The  num¬ 
ber  of  squares  was  chosen  to  cause  a  delay  of  about  100  s 
between  dot  writes.  The  results  of  such  an  experiment  per¬ 
formed  on  both  GaAs  and  Si  substrates  are  shown  in  Fig.  11. 
The  cross-sectional  profiles  are  of  the  dots  written  inside 
squares  of  various  sizes  and,  for  comparison,  isolated  dots. 
The  heights  of  the  dots  written  inside  the  squares  are  much 
less  than  those  of  the  isolated  dots.  Additionally,  as  the 
squares  are  reduced  in  size,  the  hydrocarbon  source  is  re¬ 
duced  further,  causing  the  dot  height  to  decrease.  This  is  as 
expected  from  the  hydrocarbon  surface  diffusion  picture. 

A  final  point  of  interest  concerns  the  distribution  of  con¬ 
tamination  for  small  delay  times.  In  particular,  the  contami¬ 
nation  in  and  near  the  direct  electron  beam  impingement 
region  (region  1)  is  less  than  that  of  the  surrounding  region 
(region  2)  for  both  the  line.  Fig.  5(d),  and  the  dot.  Fig.  6(d). 
This  could  arise  in  the  diffusion  picture  from  the  hydrocar¬ 
bons  being  polymerized  before  they  diffuse  into  the  center  of 
the  exposed  region.  Hence  the  contamination  would  prima¬ 


rily  be  at  the  edges  of  the  exposed  region.  For  such  a  process 
to  take  place,  the  delay  time  between  exposures  must  be  long 
enough  for  diffusion  to  occur.  Based  on  the  estimated  diffu¬ 
sion  time  of  3-10  s  determined  earlier,  this  could  be  the  case 
for  the  line  of  Fig.  5(d)  that  has  an  associated  delay  time  of 
5.2  s,  but  it  is  unlikely  that  this  could  explain  the  dot  of  Fig. 
6(d)  that  was  written  with  a  very  short  delay  time  of  22  ms. 
Another  possibility  is  that  etching  reactions  begin  to  domi¬ 
nate.  Carbon  compounds  can  be  etched  through  electron 
beam  induced  chemical  reactions  with  adsorbed  residual  gas 
molecules,  resulting  in  volatile  compounds  with  carbon.1'2 
Hydrocarbons  initially  present  in  the  directly  exposed  region 
are  polymerized  early  in  the  exposure  cycle.  As  additional 
exposures  are  made,  no  mobile  hydrocarbons  are  available 
for  further  polymerization.  The  competing  etching  reactions 
then  act  to  remove  the  contamination  previously  created. 
This  could  take  place  primarily  near  the  direct  impingement 
region  as  a  result  of  the  greater  amount  of  energy  deposited 
in  this  region.  A  final  possibility  is  that  the  exposure  causes  a 
densification  of  the  substrate  oxide.  Such  a  densification  has 
been  previously  proposed  and  observed  for  electron  beam 
exposed  silicon-oxide  surfaces.4-7  The  densification  would 
explain  the  apparent  pit  in  the  direct  impingement  region  of 
the  dot  shown  in  Fig.  6(d). 

IV.  SUMMARY 

Electron  beam  exposure  of  a  surface  can  create  a  buildup 
of  hydrocarbon  contamination  through  polymerization  of 
surface  hydrocarbon  molecules.  We  have  studied  the  topog¬ 
raphy  of  such  contamination  structures  using  atomic  force 
microscopy.  At  sufficiently  large  electron  doses,  the  struc¬ 
tures  are  found  to  be  larger  than  the  area  of  direct  electron 
beam  impingement.  Electrons  scattered  within  the  substrate 
and  then  out  through  the  surface  lead  to  the  contamination 
outside  of  the  direct  impingement  region.  The  spatial  extent 
of  contamination  beyond  the  primary  beam  area  has  been 
determined  from  cross-sectional  profiles  of  contamination 
lines  written  on  both  GaAs  and  Si  substrates.  This  extent 
compares  favorably  to  the  surface  interaction  length  pre¬ 
dicted  by  electron  scattering  theory. 

The  dynamics  of  the  contamination  growth  have  been  in¬ 
vestigated  by  writing  a  pattern  with  the  electron  beam  mul¬ 
tiple  times,  with  a  controllable  delay  time  inserted  between 
each  of  the  exposures.  For  large  delay  times,  significant  con¬ 
tamination  results  since  there  is  sufficient  time  between  ex¬ 
posures  for  the  reduced  quantity  of  hydrocarbons  in  the  ex¬ 
posed  region  immediately  following  an  exposure  to  be 
replenished  before  the  next  exposure.  The  quantity  of  con¬ 
tamination  in  this  limit  is  proportional  to  the  number  of  elec¬ 
tron  beam  exposures.  As  the  delay  time  is  reduced,  the  hy¬ 
drocarbon  surface  density  has  less  time  to  recover  and, 
consequently,  less  contamination  is  formed.  For  small  delay 
times,  the  exposed  region  can  become  completely  depleted 
of  mobile  hydrocarbons  since  there  is  insufficient  time  for 
the  transport  of  hydrocarbons  to  the  exposed  region  between 
exposures.  The  contamination  structure  in  the  direct  im¬ 
pingement  region  is  observed  not  to  grow  with  the  number  of 
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exposures  in  this  limit,  and  competing  reactions  that  reduce 
the  structure  height  may  play  a  role.  Such  dependencies  have 
been  observed  for  both  line  structures  in  which  the  delay 
time  is  dictated  by  the  line  length,  and  dots  for  which  an 
artificial  delay  time  has  been  inserted  between  exposures. 

The  observed  dependence  of  the  contamination  height  on 
the  delay  time  is  consistent  with  a  hydrocarbon  surface  dif¬ 
fusion  picture.  Within  this  picture,  the  hydrocarbons  in  the 
exposed  region  are  replenished  through  surface  diffusion  of 
hydrocarbons  from  the  surrounding  areas.  The  model  pre¬ 
dicts  the  smaller  rise  in  contamination  height  with  delay  time 
observed  with  the  line  structures  as  compared  to  that  of  the 
dots.  This  occurs  because  a  section  of  the  line  away  from  its 
ends  effectively  has  a  smaller  source  of  hydrocarbons  than 
that  of  a  similarly  written  dot.  The  line  ends  also  have  a 
greater  amount  of  contamination  than  segments  away  from 
the  ends  for  a  similar  reason. 

The  role  of  surface  diffusion  in  the  formation  of  contami¬ 
nation  has  been  further  illustrated  by  reducing  the  source  of 
hydrocarbons  to  a  region  on  the  substrate  through  electron 
beam  exposure  and  then  writing  a  dot  in  this  region.  The 
height  of  such  a  dot  is  much  less  than  that  of  an  isolated  dot, 
as  is  expected  from  the  surface  diffusion  picture. 
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APPENDIX 


In  this  Appendix  the  diffusion  equation  is  solved  for  two 
different  geometries:  one  of  rotational  symmetry  applicable 
to  the  dot  structures  and  the  other  of  translational  symmetry 
relevant  to  the  line  structures.  For  both  cases,  the  situation  is 
considered  in  which  a  single  electron  beam  exposure  reduces 
the  equilibrium  mobile  hydrocarbon  density  n0  by  a  constant 
amount  over  a  region  whose  shape  is  appropriate  for  the 
structure  of  interest.  The  diffusion  equation  is  solved  in  order 
to  determine  the  time  necessary  for  the  reduced  density  re¬ 
gion  to  reattain  the  equilibrium  value.  Consider  first  the  dot 
structure  case.  To  take  advantage  of  the  rotational  symmetry 
of  the  dot  geometry,  polar  coordinates  will  be  used, 
x  =  (r,6).  If  the  exposed  region  is  assumed  to  be  circular, 
then  n  ( r,  0 ,  t )  will  be  independent  of  6  and  the  diffusion 
equation  (4)  becomes 
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For  simplicity,  further  assume  a  region  r<d  of  reduced  hy¬ 
drocarbon  density  nr  after  the  exposure  and  no  polymeriza¬ 
tion  occurring  outside  this  region.  The  initial  condition  is 
then 


n(r,t  =  0)  =  ' 


n  o, 


0  <r<d 
r>d 


(A2) 


The  end  of  the  electron  beam  exposure  corresponds  to  t  =  0. 
The  solution  to  Eq.  (Al)  given  the  initial  condition  (A2)  is20 
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Jo 
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n(s,t=0)s  ds. 


(A3) 


where  IQ(z)  =  Jo(iz )  is  the  zero-order  modified  Bessel  func¬ 
tion.  The  density  is  minimum  at  the  center  of  the  exposed 
region  and  is  therefore  the  value  of  interest.  Setting  r  equal 
to  zero  and  using  70( 0)  =  1 ,  the  following  is  obtained  from 
Eq.  (A3): 


(A4) 


n(r  =  0,t)  =  n r  +  (nQ  —  n r)exp|  -  ■ 


The  analogous  situation  for  a  line  structure  will  now  be 
considered.  If  only  the  region  of  the  line  far  from  the  ends  is 
considered,  the  structure  has  translational  symmetry.  Using 
rectangular  coordinates,  x  =  (x,y),  will  then  simplify  the 
problem.  If  the  y  axis  is  chosen  to  be  parallel  to  the  line, 
n(x,y,t)  will  be  independent  of  y  and  the  diffusion  equation 
(4)  reduces  to  the  one-dimensional  diffusion  equation: 

dn(x,t)  d2n(x,t) 

~  (A5) 


dt 


=  D 


dx 


An  initial  condition  for  a  reduced  density  line  that  is  analo¬ 
gous  to  that  of  Eq.  (A2)  is 


n(x,t  =  0)  = 


1*1  <d 
\x\  >d 


(A6) 


The  solution  to  the  one-dimensional  diffusion  equation  (A5) 
given  the  initial  condition  (A6)  is21 


n(x,t )  = 
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/ -  I 
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=  nr+  -  (. n0-nr ) 
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where  erfc(z)  is  the  complementary  error  function 


erfc(z)  =  — 
V^r 


: 


exp(  —  s2)ds. 


(A8) 


(A9) 


Letting  x  =  0,  the  following  is  obtained  for  the  density  at  the 
center  of  the  line: 


n  (x  =  0,t )  =  n  r  +  (n  o  —  n ,.)  erfc 


(A10) 
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